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Abstract 

We present a method (the Aufbau/Abbau method) for optimizing the structure of a whole se- 
ries of clusters without making any assumptions on the structure. Subsequently, the method is 
combined with the embedded-atom method in determining the structure of the two energetically 
lowest isomers of Ni^r clusters with N up to 150. Finally, various analytical descriptors are intro- 
duced that are used in studying the overall shape of the clusters, their structure and stability, and 
possible growth and dissociation processes. 
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About two decades ago the experimental observation and theoretical description of so- 
called magic numbers in clusters of simple metals (see 313]) provided an excellent example 
of how experiment and theory could supplement each other in the development of the un- 
derstanding of materials properties. The theoretical studies for any but the smallest clusters 
were based on the jellium model, i.e., the internal structure of the cluster was completely 
neglected and the core electrons and the nuclei were smeared out to a homogeneous jellium 
in whose field the valence electrons moved. Moreover, in the first studies the jellium was 
supposed to be spherical and to have sharp boundaries, although these assumptions in later 
studies have been relaxed partially (see ^|). 

For theoretical studies of any other metal clusters M^r than those of simple metals, one 
encounters the serious problem that when the valence electrons are anything but completely 
delocalized, the (spherical) jellium model may not provide a realistic approximation. This 
is, e.g., the case for transition-metal clusters where the nearest-neighbor interactions be- 
tween the atoms are important, so that the precise structure of the cluster as well as the 
inhomogeneities of the electron density are to be included in any realistic description. Then, 
furthermore, the fact that the number of structural degrees of freedom (3A^ — 6 or 3A^ — 5) 
together with the fact that the number of metastable structures grows very fast with 
makes detailed systematic studies of how the properties depend on for any but the abso- 
lutely smallest clusters extremely demanding. On the other hand, the further development 
of this field requires a precise understanding of how the properties of interest depend on 
the structure, size, and composition of the clusters, which, unfortunately, is very difficult to 
achieve theoretically or experimentally. 

Here we shall report results of a theoretical study of the structural and energetic proper- 
ties of NIat clusters for all A^ up to 150. Our calculations have made no assumptions about 
the structure of the clusters whereby they differ from most other studies of these systems 
(see. e.g.. fflffl and references therein) of which the previously most detailed and un- 
biased study^ considered only up to A^ = 80. Since we have had to perform often several 
1000 structure-relaxations for a given A^ without making any assumption on structure or 
symmetry, parameter-free methods cannot be applied, and, therefore, we have employed the 
embedded-atom method (EAM) in the calculation of the total energy for a given structure. 
Moreover, we have developed an Aufbau/Abbau method for the structure optimization. This 
method is partly inspired by experimental conditions where clusters often are produced by 
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adding atoms individually to a seed, although the experiments often take place at conditions 
far from equilibrium and from T = K. A preliminary description of our approach has been 
presented elsewhere,—*^ but compared with that work the present one is much more accurate 
(through a much more detailed structure optimization that also includes more isomers), and, 
in addition, we shall here present various theoretical tools with which we can extract the 
fundamental structural and energetic properties of these clusters. 

In our Aufbau/Abbau method we consider clusters with L and L + P atoms with P ~ 
5 — 10. For each of those we study a set of randomly generated structures. Using a quasi- 
Newton method relaxed local-total-energy minima are identified and the structures of the 
lowest total energy selected. One by one atoms are added to the structure with L atoms 
(many hundred times for each size), and the structures are relaxed. Parallel, one by one 
atoms are removed from the structure with L + P atoms. From this two series of structures 
for L < N < L + P those structures of the lowest energies are chosen and these are used 
as seeds for a new set of calculations. First, when no lower total energies are found, it is 
assumed that the structures of the global-total-energy minima have been identified, and we 
proceed to larger clusters. In some few cases, extra high-symmetry structures were added 
as possible structures, which in one case (A^ = 75) resulted in an otherwise unidentified 
lowest-energy structure. According to our experience, for a given A^ we have to consider 
about 500 randomly generated structures for the smallest values of A^ and about 1000 for 
the largest values. Moreover, about 3 Aufbau/Abbau cycles are needed for the smaller A^ 
and up to 10 for the larger A^. 

In comparison with our earlier study, the method has been extended so that by keep- 
ing track of more structures of the lowest total energies, we obtain the energetically lowest 
isomers for each size, and a much more detailed search has led to the identification of new 
structurals for several A^. Since our approach (as any other unbiased method for structure 
optimization of larger, low-symmetric systems) requires a large number of total-energy cal- 
culations for different, larger structures, it cannot be combined with accurate parameter-free 
electronic-structure methods, but has to be based on approximate schemes for calculating 
the total energy. Here, we combine the Aufbau/Abbau method with the EAM method for 
the latter purpose. 

The basic idea of the EAMiSiiiJ^ii^ is to consider every atom as an impurity embedded 
in a host provided by all other atoms. The tails of the electron densities of the host atoms 
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inside the guest atom lead to a modification in the energy of this atom, Fi{p^), where i labels 
the atom, and is the electron density from the host atoms. Usually this is evaluated at 
the site of the nucleus of the ith atom. The remaining part of the binding energy of the 
system of interest is approximated through pair potentials, \J2j^i4'ij{Rij) with Rij being 
the distance between atoms i and j. The embedding energies Fi and the pair potentials 
(pij have been obtained by fitting to theoretical and experimental information on different 
infinite systems.— The EAM is conceptually simple, but it grasps the major parts of the 
interatomic interactions and gives accurate structural and energetic information on various 
extended and truncated systems.— 

Before discussing our results in detail we shall address the question of the accuracy of 
the EAM method. As such, it is constructed first of all for describing cohesion in extended 
metallic systems like crystals without or with defects and surfaces. Due to the complexity of 
the present study, it is not possible to apply a parameter-free method, and simpler methods 
like the EAM have to be used, but nevertheless it may be asked whether the EAM can be used 
for the, after all, relatively small systems of our present study. We shall therefore compare 
our calculated structural properties for nickel clusters with = 2, 3, 4, 5, 6, 7, 8, and 
13 with those of accurate electronic-structure calculations using different density-functional 
methodsi^^*i^*i^*i2»i^ Table |l] list some of the key quantities in the structural characterization 
of the different clusters as obtained with the different methods, and in Fig. ^ our optimized 
structures are shown. It is seen in Table n that our results agree well with those of the other 
studies, both concerning the structure itself and concerning the interatomic distances. Most 
of the discrepancies are related to smaller distortions that may be due to electronic effects 
(e.g., Jahn- Teller distortions) that are not included in the EAM method. Thus, we conclude 
that our approach is accurate. 

In Fig. 121 we show our calculated binding energy per atom for nickel clusters with N 
ranging from 1 to 150. As is generally found, the binding energy for clusters is first of all 
a monotonously increasing function of with, however, some extra structure. In order to 
identify the latter more clearly, we show in the figure also the stability function S{N) = 
Etot{N + 1) + Etot{N — 1) — 2Etot{N) which has maxima for particularly stable structures. 
These are also marked in the figure. As we shall see, some of the particularly stable structures 
correspond to particularly symmetric ones. Particularly stable structures are found for A^ 
= 13 (the most pronounced one), 19, 23, 39, 55, and 147. 
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There exists some previous experimental and theoretical studies of the properties of Nijv 
clusters, of which many focus on the inner structural aspects (including symmetry, see, e.g., 



22ll23lJ24l ]). The most detailed study of the energetics of Ni^v clusters is the theoretical 
one of Doye and Wales'^ who used the Sutton-Chen potential in calculating structural and 
energetic properties of Ni^v clusters with up to 80. They find particular stability for = 
13, 38, 50, 55, 64, 71, 75, and 79, which is in partial agreement with our results. As discussed 
by Parks et alM^ the structures given by Doye and Wales seem to deviate somewhat from 
experimental findings, suggesting that the Sutton-Chen potential is not sufficiently accurate 
for such studies. In a tight-binding molecular-dynamics study, Luo studied clusters with 
up to 55 atoms. He found most stable structures for A^ = 21, 31, 35, 38, 40, 47, 50, and 
55, only partly in agreement with the results of Doye and Wales and of us. On the other 
hand, Montejano-Carrizales et al^ used the EAM in analysing the so-called umbrella growth 
models for icosahedral Ni clusters. Although they thereby make certain assumptions on the 
structure their A^ values for particularly stable structures are in general in good agreement 
with ours and with those of Doye and Wales. This leads support to the umbrella growth 
model. This conclusion is further supported below. 

The special stability of the cluster sizes mentioned above is supported by comparing the 
energies of the two energetically lowest isomers. First for A^ = 6 and onwards we identify 
two different isomers, and their total-energy difference is shown in the lowest panel in Fig. 
12] Here, the most significant features correlate closely with those of the stability function. 

Besides the stability a number of other issues have attracted significant interest in cluster 
physics. These include the structure of the individual clusters, whether growth modes can 
be identified, and how the clusters will dissociate. We shall now address these issues. 

For a given cluster NIat we calculate the eigenvalues laa of the matrix containing the 
moments of inertia Yu^=i Siti with s and t being the x, y, and z in a coordinate system with 
the origin at the center of mass of the cluster. For a spherical jellium laa oc N^^^, and 
therefore we show (Fig. IH} laa/N^^^- If the three eigenvalues are identical the overall shape 
of the cluster is approximately spherical. Such clusters are marked with points in the lowest 
row in the uppermost part of the figure and are found for A^ = 1, 4, 6, 13, 26, 28, 38, and 
55, which partially correspond to particularly stable clusters; cf. Fig. |21 When two of the 
eigenvalues are small and only one is large, the cluster is essentially prolate (cigar-shaped), 
whereas it with two large eigenvalues and one small one is oblate (lens-shaped). In the figure 
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we have marked which of the clusters have these two types of overall shapes. Except for 
the absolutely smallest clusters, we have that, if the cluster structures can be considered as 
obtained by adding atom by atom to a largely fixed core, the clusters will either be close 
to spherical (so that several changes between prolate and oblate structures are possible), or 
one or the other of the two structure types will dominate in large A^-intervals. However, 
the results of Fig. El show that except for above about 137, the eigenvalues laa are very 
different, so that the clusters have pronounced prolate or oblate shapes. Moreover, several 
changes between these two structure types occur in the range of we have considered. 

A more detailed way of analysing whether the cluster with A^ atoms can be considered 
as formed by adding a single atom to the cluster with A^ — 1 atoms is provided by defining 
similarity indices as follows. We calculate the (A^ — 1)(A^ — 2)/2 interatomic distances dP^^ (in 
a.u.) for the (A^ — l)-atom cluster and sort them. Subsequently, we consider all the possible 
(A^ — l)-atomic parts of the A^-atom cluster, and for those we also calculate and sort the 
interatomic distances dij. The smallest value of qi with (^~^H^~^) g2 _ ^^"^^(rfO. — dijY 
defines then the similarity index, si = (1 + Q'l)^^, which is 1 for a cluster that is obtained 
by adding one atom to the previous one, and approaches for very different structures. 
Another similarity index S2 is obtained by instead of using the interatomic distances we use 
the distances of the individual atoms to the center of mass of the cluster. Denoting these 
so-called radial distances r° (in a.u.) for the (A^ — l)-atom cluster and (in a.u.) for the 
(A^ — l)-atom fragment of the A^-atom cluster, we consider the smallest value of q2 with 
{N - l)ql = E.=T'(r° - nf, and define = (1 + ^2)"^ 

These two functions are shown in Figs. Efa) and (b). Up to A^ around 40 both functions 
possess a number of peaks indicating that no simple growth mechanism can be identified. 
This is in accord with the experimental observations of Pellarin et ali^ who moreover suggest 
that above that the umbrella growth model agrees with their observations, in agreement with 
the results of Montejano-Carrizales et al^ However, both of our similarity indices find for 
larger A^ peaks at A^ = 39, 51, 63, 73, 80, 87, 105, and 138, implying that also here size- 
dependent structural differences occur. Thus, in total we find an overall support for the 
umbrella growth model although for many values of A^ this model does not give an exact 
description of the structure. 

In order to obtain a quantitative estimate of the similarity of the structure for a given A^ 
with that of a fragment of the crystal we proceed as follows. We consider a large, but finite 
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fragment of the fee erystal strueture and plaee a eenter (arbitrarily) at (0,0,0), (|a,0,0), and 
(|a,|a,0) (with a being the length of the unit cell). The set of distances from the atoms to 
any of these centers, is sorted and compared with the radial distances for a given A^- 

atom cluster, Nqj^^ = YliLi{rf''^ — riy. This defines the fee similarity index, Sfcc = (1+gfcc)^^- 
A related icosahedral similarity index Sico is defined by using the radial distances for the 
(relaxed) = 309 icosahedral cluster in comparing with the different clusters. The two 
functions are shown in Figs, ^c) and (d). It is seen that in the range 73-79 the structure is 
related both to that of a fragment of the fee crystal with, in addition, the center at (0,0,0) 
and to the icosahedral structure. This indicates that in exactly this size range more different 
structural motifs are present. Two further peaks can be identified in Fig. EJ^c): at = 28 
which corresponds to the center at (|a,0,0) and at = 38 which corresponds to the center 
at (^a,|a,0). The lower panel shows that the icosahedral clusters are important around 
A^ = 13, between A^ = 51 and A^ = 55, and above A^ = 138. We finally notice that the 
structural characteristics often are not isolated to singular values of A^ but can be identified 
also around the peaks. 

In order to get insight into possible fragmentation channels we analyse the calculated 
energies i^tot as a function of A^. Neglecting kinetic effects (this is a crude approximation 
and may have severe consequences for our conclusions), we seek for a given A^ that value of 
K for which the dissociation energy Etot{N) — [Etot{K) + Etot{N — K)] is lowest. It turned 
out that, without exception, this value of K was always 0, indicating that any cluster is 
stabler than its dissociation products. This result is a consequence of the fact that Etot{N) 
is upward curved, as indirectly can be seen in the upper part of Fig. |21 

Neglecting the K = solution we found the results shown in Fig. |S| Except for one 
case, K was found to be either 1 or 2, indicating that the clusters upon fragmentation will 
dissociate in a highly asymmetric way by splitting off a single or a pair of atoms. The 
exception is that of A^ = 16 for which K = 3 was found as a consequence of the particular 
stability of the A^ = 13 cluster. The dissociation energies, also shown in Fig. El are seen to 
lie around 4 eV for K = 1 and 2 and around 5 eV for if = 3 for the larger clusters. 

Concluding we have presented results of a theoretical study of structural and energetic 
properties of Ni^r clusters with all A^ up to A^ = 150. We stress that this study is the first 
where a systematic study of energetics, growth, dissociation, and structure of not only the 
absolutely smallest clusters has been performed without more or less severe assumptions 
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of the structure. Our structure optimization was performed using our own Aufbau/Abbau 
method, and due to the large number of total-energy calculations that was necessary we 
used the approximate EAM for this purpose. As the main result we obtained the structure 
and the total energy as functions of A^, and these results were then subsequently used 
in identifying particularly stable clusters, in classifying the overall shape of the clusters, 
in discussing whether a simple growth mechanism could be applied (this was found not 
to be the case), and in discussing possible fragmentation channels. In order to extract 
useful information from the obtained structural and energetic properties we introduced a 
number of quantitative measures and subsequently applied them on our results. For the 
first time it was shown that for 73 < A^ < 79 the clusters possess more different types of 
structural characteristics, and, moreover, our theoretical study is the first such to show that 
for A^ = 147 a magic number occurs for a multilayer icosahedron structure. This result is in 
perfect agreement with mass-spectroscopical experiments of Pellarin et alM 
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Dooh 

2.12 

Dzh 
2.25 

Td 
2.32 

2.35 

Oh 
2.36 

2.39 

2.38 



Dooh 

2.17 



Dooh 

2.13 
2.26 
2.33 
2.36 

C^ (~ Oh) 

2.40 



2.32 



C2v (~ -Ds/i) 
2.28 

2.28 



2.37 



Dooh 

1.99 

2.15 

D2d I 
2.17 / 2.10 

2.25 

Ah 
2.33 



Oh 
2.16 

2.23/2.34 



2.36/2.48 2.41/2.53 

TABLE I: Our optimized structures in comparison with those of ab initio density-functional cal- 
culations. In two cases {N = 5, Ref. Ilal a nd N = 7, Ref. ,17]) we compare with the second-lowest 
structure, and in another (A^ = 4, Ref. 19|) we give their results for the two energetically close min- 



ima. The second row gives averaged bond lengths in A. With we mark 



The experimental bond length for the dimer is 2.15—2.20 A (Refs 



approximate symmetries. 
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FIG. 2: The binding energy per atom for the energeticahy two lowest isomers (upper part), the 
stabihty function (middle part), and the total-energy difference for the two energetically lowest 
isomers (lowest part) for Nijv clusters as functions of N. 
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FIG. 3: The normalized eigenvalues of the matrix containing the moments of inertia as functions 
of N. The three solid curves gives the eigenvalues, and the dashed curve their average. In the 
upper part of the figure a simple estimate of the overall shape of the clusters is given: either 
being spherical (marked with points in the lowest row), cigar-shaped (points in the middle row), 
or lens-shaped (points in the uppermost row). 
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FIG. 5: The fragmentation channels for the Nijv clusters as functions of N. The upper panel 
shows the most probable size of the smallest split-off cluster, whereas the lower panel shows the 
dissociation energy for splitting off 1 (lowest curve), 2 (middle curve, shifted upwards by 2 eV), or 
3 atoms (upper curve, shifted upwards by 4 eV). 
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